Optically detecting the quantization of collective atomic motion 
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We directly measure the quantized collective motion of a gas of thousands of ultracold atoms, 
coupled to light in a high-finesse optical cavity. We detect strong asymmetries, as high as 3:1, 
in the intensity of light scattered into low- and high-energy motional sidebands. Owing to high 
cavity-atom cooper at ivity, the optical output of the cavity contains a spectroscopic record of the 
energy exchanged between light and motion, directly quantifying the heat deposited by a quantum 
position measurement's backaction. Such backaction selectively causes the phonon occupation of 
the observed collective modes to increase with the measurement rate. These results, in addition 
to providing a method for calibrating the motion of low-occupation mechanical systems, offer new 
possibilities for investigating collective modes of degenerate gases and for diagnosing optomechanical 
measurement backaction. 

PACS numbers: 07.60.Ly,42.50.-p,42.65.Dr 



Quantized motion leads to a large asymmetry in the 
spectrum of light scattered by a ground-state oscillator. 
Such asymmetry is most commonly observed from micro- 
scopic oscillators, such as electrons bound within atoms 
and molecules, single neutral atoms [1], or few-ion ensem- 
bles [2-4]. In contrast, quantum aspects of the motion of 
massive, many-atom oscillators are typically obscured by 
thermal noise and high phonon occupation. These ob- 
jects, therefore, usually modulate the spectrum of light 
in a classical, symmetric manner. 

Controlling and measuring the motion of macroscopic 
objects at levels sensitive to quantum efTects will be crit- 
ical for operating gravitational- wave detectors [5] , verify- 
ing the correspondence principle at macroscopic scales 
[6, 7], and realizing protocols that mechanically store 
and exchange quantum information [8, 9]. Such goals 
are being pursued actively using cavity optomechanical 
systems, wherein the motion of an object with mass rang- 
ing from attograms to kilograms is observed via its cou- 
pling to an electromagnetic cavity [10]. Implementations 
of cavity optomechanics using the collective motion of 
atomic gases [11, 12] have demonstrated sensitivity to 
quantum optical force fluctuations, leading to the obser- 
vation of optical squeezing from ponderomotive interac- 
tions [13] and backaction from a quantum-limited posi- 
tion measurement [14]. 

In this work, we use a high-finesse optical cavity to 
detect the coherent, asymmetric scattering of light by 
collective modes of motion of a trapped atomic gas, oc- 
cupied with as few as 0.5 phonons. Observations have 
previously been made of optical emission asymmetries 
from individual atoms [1] and ensembles of up to 14 ions 
[2-4] , and of asymmetric absorption by a nanomechanical 
solid-state resonator [15]. In this work, we measure the 
coherent scattering of light from the collective motion of 
thousands of ultracold atoms. The scattering asymme- 
try acts as a self- calibrating thermometer for the atoms' 



collective phonon occupation. Moreover, owing to high 
cavity-atom cooperativity and thermal isolation in our 
system, the cavity mode acts as the dominant channel 
for energy fiux to our mechanical system. The spec- 
trum of light emitted from the cavity therefore serves 
as a record of the energy exchanged between motion and 
the light field. We demonstrate that this energy trans- 
fer represents the necessary minimum diffusive heating of 
a continuous backaction-limited quantum position mea- 
surement. 

Our experiment (Fig. 1) begins with an ensemble of 
4,000 ultracold ^^Rb atoms. The atoms are trapped 
in a few adjacent minima of a one-dimensional optical 
standing-wave potential, formed by 850-nm-light light 
resonating within a high-finesse Fabry-Perot optical cav- 
ity, and detuned far from atomic resonance. The cur- 
vature at each potential minimum corresponds to an 
oscillation frequency, along the cavity axis, of Um = 
27r X 110 kHz. 

We create an optomechanical coupling, which is lin- 
early sensitive to atomic position, by trapping the atoms 
at locations with strong intensity gradients of 780-nm- 
wavelength probe light, which also resonates within 
the cavity [16]. The probe light is detuned from the 
atomic D2 transition by an amount Aca, ranging from 
— 27r X 12 to 70 GHz. At such detunings, the atomic gas 
acts as a position-dependent refractive medium, lead- 
ing to an interaction Hamiltonian Hint = ^i^^i/^HO, 
where h is the probe's photon number operator and Zi is 
the position operator of atom i [14]. Here gi represents 
the change in the cavity's resonance frequency as atom i 
is displaced by one single- atom harmonic oscillator length 
^HO, equivalent to ^jTil2muJra for rubidium mass m. 

The above may be rewritten [11, 12] as a collective 
interaction described by ^int = QomnZ / Z^o^ where 
Z = ^iQiZi/ ^iQi is a collective position operator, 
Z^io = ^Ro/V^ is the collective oscillator length for 
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FIG. 1: Schematic of the experiment. (A) An ensemble of 4,000 ultracold ^^Rb atoms is trapped in an optical standing wave 
within a high-finesse optical cavity. The ensemble's center-of-mass motion coherently scatters light from cavity probe light, 
detected using a balanced heterodyne receiver. (B) Probe spectrum in the cavity. A strong resonant tone (green) applies 
backaction to and acquires sidebands from the collective atomic motion. A weak detuned tone (gray) is used for locking the 
probe frequency with respect to the cavity, and does not significantly affect the oscillator, neither via incoherent backaction 
nor via dynamical cooling. (C) A harmonic oscillator in its ground state can only extract energy from the optical field, leading 
to an asymmetry in the resonant probe's Raman sidebands, indicative of the mean phonon occupation number v. 



Na atoms, and ^om = is the collective op- 

tomechanical coupling rate. In our system, the gi are all 
approximately equal, and Z is nearly equivalent to the 
center-of-mass position ^com = ^i^i/^a- In this case, 
^om = V^^o ^^Ho/ Aca, with siuglc-atom cavity-QED 
coupling rate go (27rxl2.4 MHz, averaged over the intra- 
cavity atomic distribution) and probe wavenumber k. 

Unlike solid-state resonators, which are coupled to 
their material environments, the atoms here are isolated 
in vacuo, and the thermal "bath" for the measured collec- 
tive oscillation consists of the remaining collective modes 
of atomic motion. After raising the trap potential to 
its final strength, the bath has an initial temperature 
of 2.9 /iK, corresponding to a mode occupation of 0.54 
phonons. 

The atomic motion scatters cavity photons via col- 
lective virtual electronic excitations, causing transitions 
between collective motional energy levels {i.e., collec- 
tive Raman scattering). In this process, the collective 
atomic oscillator can absorb a mechanical energy quan- 
tum, thereby downconverting the emitted photon by uOm 
(Stokes scattering). If the oscillator is not in its ground 
state, it can alternatively upconvert the photon by cj^, 
thereby losing an energy quantum (anti- Stokes scatter- 
ing). By Fermi's golden rule, the rates for upconver- 
sion and downconversion, from an oscillator in its uth 
excited state, are proportional to \{u\Z\u—l)\'^ ex u and 
I P (X u -\- 1, respectively. A ground-state os- 
cillator {u = 0) will thus scatter only to the Stokes 



sideband, with more symmetric scattering indicative of 
higher phonon occupation. 

The Fourier spectrum n{uj) of photons emitted from 
the cavity, at a frequency u relative to the probe, and 
normalized by the discrete Fourier window time, can be 
calculated by solving a set of quantum Langevin equa- 
tions [15, 17, 18], yielding 
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where u is the oscillator's mean phonon occupation. Here 
we have parameterized the photon scattering rate by the 
dimensionless optomechanical cooperativity [19] Com = 
4n^Qj^/r^/^, which combines ^om, the mechanical damp- 
ing rate F^, the mean intracavity probe photon number 
n, and the cavity half-linewidth (27r x 1.82 MHz). In 
our atoms-based optomechanical system. Com n/A^^, 
and the cooperativity can be tuned over many orders of 
magnitude by varying the probe intensity and Aca- 

We measure the Stokes asymmetry by integrating the 
optical power P± scattered to frequencies near icj^, 
detected using a balanced heterodyne receiver. The 
detector measures a power spectral density 6'het(^) = 
^SN [1 + e{no{uj) + n(co'))/2], where 5'sn is the mean shot- 
noise spectrum as measured by the detector and no(co') 
is the spectrum of technical noise, which contributes 
less than 1% to the observed spectrum at icj^- The 
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FIG. 2: Asymmetric optical scattering from quantum collec- 
tive motion. (A) Photon spectrum of light exiting the cav- 
ity, scattered by collective atomic motion. Shown are the 
Stokes sidebands (left panels, red) and anti-Stokes sidebands 
(right panels, blue) at various cooperativities (from top to 
bottom. Com — 9.6, 1.9, 0.4), together with the prediction of 
Eqn. (1), plotted vs. cyclic frequency f = uj/2ti. The mechan- 
ical linewidths and resonance frequencies are fit to the data, 
but the peak heights are not free parameters. (B) Measured 
phonon occupation vs. cooperativity. The collective mode 
occupation (green circles) increases by Com/2, according to 
the zero-free-parameter measurement backaction theory (red 
line, shaded region indicates 68% systematic confidence in- 
terval). In contrast, the r.m.s. single-atom axial occupation, 
measured using time-of-flight thermometry of the gas (yel- 
low diamonds, measured at 1 ms of probing), remains near 
its initial value during the measurement (gray line indicates 
theoretical prediction) . 



quantum efficiency for measuring intracavity photons is 
e — 16% ±2%. For each run of the experiment we mea- 
sure the spectral density Sa with atoms and 5*0 without 
atoms (which measures no), and then measure 5'sn by 
extinguishing the probe beam. The unitless photon spec- 
trum is then n(co') = 2(8 a — So)/eSsN- In order to avoid 
excessive accumulated heating of the atomic gas, we ac- 
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FIG. 3: Energy exchange between light and motion. (A) 
Energy-exchange spectral densities for various probe coopera- 
tivities (from top to bottom. Com = 9.6, 3.2, 1.6). Traces are 
vertically offset for clarity (by 12, 2, x 10"^'^ W/Hz). (B) 
Heating rate due to the resonant probe (magenta squares), 
the weak locking tone (gray circles), and the summed heat- 
ing from both tones (green diamonds), together with zero- 
free-parameter measurement-backaction theory (gray solid 
line), plotted as a function of the collective coupling rate 
= FmCom. (C) Temperature of the atomic gas vs. 
the product of the collective coupling rate and the total inter- 
rogation time (yellow diamonds at 5 ms, cyan circles at 9 ms). 
Also shown is the theoretical increase in the bath temperature 
(gray hne). 



quire data for only 5 ms each run of the experiment, after 
which we observe the mechanical resonance to broaden 
significantly. 

Measured spectra are shown in Fig. 2. At the low- 
est optomechanical cooperativites {viz., a weak probe de- 
tuned far from atomic resonance), the probe contributes 
only minimally to the phonon occupation. Rather, the 
collective mode occupation should approach the average 
thermal occupation (0.54 ± 0.02 phonons, measured via 
time-of-flight thermometry). We measure a large Stokes 
asymmetry: P_/P+ = 3.0 ± 0.8, corresponding to a 
phonon occupation u = P^/{P_-P^) = 0.49±0.10. We 
note that such a "Stokes thermometer" is self-calibrating, 
in that no experimental parameters are involved in ex- 
tracting the phonon occupation number. 

The mechanical resonance widths of the observed spec- 
tra are generally larger than the true mechanical damping 
rate, which we find by probing to the blue of cavity res- 
onance and finding the phonon lasing threshold [18, 20], 
at which point the mechanical damping rate is equal to 
the optomechanical amplification rate. In our system, we 
find a damping rate of 27r x 1.5 kHz. Comparing to the 
observed inhomogenous linewidth of 27r x 3 to 4 kHz indi- 
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cates that Z dephases into approximately 2 to 3 collective 
normal modes, each containing 1,500 to 2,000 atoms. We 
suspect that this dephasing is due mostly to slight dif- 
ferences in trap curvature in adjacent potential minima 
arising from the probe beam [16], and partly to the an- 
harmonicity of the standing-wave optical trap. 

We next apply the sideband thermometer to probe the 
effects of measurement backaction. A quantum mea- 
surement of position must be accompanied by a cor- 
responding motional disturbance [21]. In this experi- 
ment, where cj^ ^ we expect an increase in the 
phonon occupation by an amount Com/2. As we in- 
crease the cooperativity, the sideband asymmetry de- 
creases as theoretically predicted. Finally, we can com- 
pare the collective mode's mean phonon occupation 
^ = ((Zl^^i^^/^om)^) /2^|o - V2 to an upper bound 
on the mean single-atom phonon occupation Va = 
Si ^^T) /^^a-s^HO ~ V25 extracted by measuring the gas 
temperature using time- of- flight thermometry. The up- 
per bound is derived by assigning the gas's entire tem- 
perature increase after 1 ms of probing (equivalent to 
many motional equilibration times) to the mean single- 
atom motion along the cavity axis. The large discrepancy 
observed between the optically detected and single-atom 
occupations (Fig. 2b) highlights the fact that our detec- 
tor senses the collective motion of the gas, rather than 
the motion of individual atoms. 

The optical spectrum moreover serves as a record of 
the energy exchanged between light and motion. A pho- 
ton recorded at frequency uj indicates the emission of 
an energy Tiu from the atoms into the cavity fleld. The 
spectral density of energy-exchange is ?iujn{uj) (units of 
W/Hz), and the total power passing from the light to 
motion is [14, 17] 



Pom = ^ Tiuj n{uj)duj : 



-L m 



Heat-exchange spectra for several values of the cooper- 
ativity, as well as heating rates taken from integrating 
over the sidebands, are shown in Fig. 3. The heating 
rates agree well with the prediction of measurement back- 
action. We additionally correlate the total energy ex- 
changed with the atoms to the increase in the (flnite) 
bath's temperature. The bath temperature should in- 
crease by an amount given by equating the backaction 
heat with the temperature integral of the gas's heat ca- 
pacity, which we calculate using Bose-Einstein statistics 
for axial motion and the ideal gas law for radial motion. 
We flnd quantitative agreement to this theory as we vary 
both the cooperativity and the probe duration. 

In this work, we have demonstrated the quantization 
of the collective motion of thousands of atoms, observing 
Stokes asymmetry and zero-point motion. The Stokes 
asymmetry provides a self- calibrating thermometer for 
low-occupation collective modes. We have in addition 
observed the spectrum of energy exchanged between light 



and collective atomic motion, spectroscopically identify- 
ing backaction heating from a quantum position measure- 
ment. While our system measures center-of-mass motion, 
other modes [12, 15, 16] of physical interest could be ad- 
dressed by tailoring the light-motion interaction. For ex- 
ample, quadrupole [22-24] or scissors [25] modes could be 
sensed using quadratic optomechanical coupling [16], al- 
lowing for precise measurements of the effects of interac- 
tions, superfluidity, and viscosity in degenerate Bose and 
Fermi gases. Sideband spectroscopy of phonon modes in 
spatially extended gases [12] provides the means to study 
thermodynamics in static and driven systems [26]. 

The authors acknowledge support from the AFOSR 
and the NSF. 



* Electronic address: nbrahms@berkeley.edu 

[I] P. S. Jessen et a/., Phys. Rev. Lett., 69, 49 (1992). 

[2] F. Diedrich, J. C. Bergquist, W. M. Itano, and D. J. 

Wineland, Phys. Rev. Lett., 62, 403 (1989). 
[3] R. Islam et al, Nat. Commun., 2, 377 (2011). 
[4] T. Monz et al, Phys. Rev. Lett., 106, 130506 (2011). 
[5] K. Goda et al, Nat. Phys., 4, 472 (2008). 
[6] W. Marshall, C. Simon, R. Penrose, and D. Bouwmeester, 

Phys. Rev. Lett., 91, 130401 (2003). 
[7] O. Romero-Isart et al, Phys. Rev. Lett., 107, 020405 

(2011). 

[8] K. Stannigel, P. Rabl, A. S. S0rensen, P. Zoller, and M. D. 

Lukin, Phys. Rev. Lett., 105, 220501 (2010). 
[9] Q. Lin et al, Nat. Photon., 4, 236 (2010). 
[10] T. J. Kippenberg and K. J. Vahala, Science, 321, 1172 

(2008). 

[II] S. Gupta, K. L. Moore, K. W. Murch, and D. M. 
Stamper-Kurn, Phys. Rev. Lett., 99, 213601 (2007). 

[12] F. Brennecke, S. Ritter, T. Donner, and T. Esslinger, 
Science, 322, 235 (2008). 

[13] D. W. C. Brooks et al, "Ponderomotive light squeez- 
ing with atomic cavity optomechanics," (2011), arXiv: 
quant -ph/ 11 07. 5609 . 

[14] K. Murch, K. Moore, S. Gupta, and D. Stamper-Kurn, 
Nat. Phys., 4, 561 (2008). 

[15] A. H. Safavi-Naeini et aL, "Measurement of the quan- 
tum zero-point motion of a nanomechanical resonator," 
(2011), arXiv:quant-ph/1108.4680 . 

[16] T. Purdy et al, Phys. Rev. Lett., 105, 133602 (2010). 

[17] F. Marquardt, J. P. Chen, A. A. Clerk, and S. M. Girvin, 
Phys. Rev. Lett., 99, 093902 (2007). 

[18] See supplemental material. 

[19] J. D. Teufel et aL, Nature, 475, 1 (2011). 

[20] K. Vahala et al, Nat. Phys., 5, 682 (2009). 

[21] A. A. Clerk, S. M. Girvin, F. Marquardt, and R. J. 
Schoelkopf, Rev. Mod. Phys., 82, 1155 (2010). 

[22] S. Stringari, Phys. Rev. Lett., 77, 2360 (1996). 

[23] A. Altmeyer et a/., Phys. Rev. Lett., 98, 040401 (2007). 

[24] C. Cao et a/.. Science, 331, 58 (2011). 

[25] O. Marago et al, Phys. Rev. Lett., 84, 2056 (2000). 

[26] K. Baumann, C. Guerlin, F. Brennecke, and 
T. Esslinger, Nature, 464, 1301 (2010). 



5 




Frequency (kHz) 



FIG. SI: Analytic evaluation of the photon spectrum when 
probing at Ape = k. Two cases, each having the same total 
optomechanical cooper at ivity, are considered: A system with 
two normal modes (left), each damped at Fm = 27rx2 kHz and 
also dephasing at a rate 27r x 2 kHz, and a system with one nor- 
mal mode (right), damped at Fm = 27rx4 kHz. Each are plot- 
ted for increasing photon number n = 0.01, 0.03, 0.1, 0.3, 0.55, 
with Qom = 27r\/Na X 3 kHz. For each, the spring shift gives 
the total optomechanical coupling, and the phonon lasing 
threshold gives the homogenous damping rate. 



SUPPLEMENTAL MATERIAL 

Photon spectrum: For a single measured normal 
mode, an optomechanical system can be described by 
the Hamiltonian 



./out • 



(SI) 



Here c^cav is the cavity resonance frequency, d is the cav- 
ity field operator, and b is the mechanical field operator. 

We solve the system's dynamics in the Fourier do- 
main, writing 0{uj) to represent the Fourier transform - 
with Fourier bandwidth ujbw ~ of an operator 0{t). We 
find the spectrum n{uj) = cjBw(^lut(^)^out(^)) of pho- 
tons exiting the cavity by solving the system's quantum 
Langevin equations in the Heisenberg picture [1-3]: 



cal bath. Using the cavity output boundary condition 
^in + ^out = V^d, assuming r^<Cco'^, and applying 
Parseval's theorem {Yl b^booBW = b^bAt, for sampling 
interval At) yields Eqn. (1) of the main article. 

Damping rate: In order to measure the native damp- 
ing rate of the collective atomic motion, we measure the 
system's phonon lasing threshold. This is done by prob- 
ing the cavity a detuning Ape = above the cavity's opti- 
cal resonance. Under this condition, the optomechanical 
spring [1, 3] shifts each normal mechanical mode's reso- 
nance frequency to a new value ujg- In the limit uj^n <^ 



i + 4ngl^A,J{^^ + Aj,). (S4) 

The spring shift is also accompanied by an optomechani- 
cal amplification rate arising from cavity feedback, equal 
to 



2k 



amp 



A2 



Sng[ 



(«2 + A2 )2 



(S5) 



When this amplification rate equals the native damping 
rate of the system, the system crosses a "phonon lasing" 
threshold; for Famp > Tm? the oscillation becomes unsta- 
ble and grows dramatically. F^ can therefore be found 
by tuning ujs until the system crosses the lasing thresh- 
old. Fig. SI shows, by solving the Langevin equations for 
one- and two-mode mechanical systems, how this tech- 
nique can be used to differentiate between a system with 
a single measured normal mode and a system with mul- 
tiple measured normal modes. 



-iujd{uj) 
-iujb{uj) 



2k: ain(cj) — K:d{uj) 
igomVn{ib\-uj)^b{uj)), (S2) 

Tm ^in 

igoinV^{d\-uj) + a(cj)), (S3) 



where d-^^iuj) and aout(^) are the cavity input and output 
fields, respectively, and 6in(co') represents the mechani- 
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